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ABSTRACT 


This bulletin is the third in a series of three bulletins 
dealing with a number of oscillatory problems. In the 
first two bulletins emphasis was placed on steady-state 
oscillations, while in this bulletin the discussion is con- 
cerned with nonsteady-state or transient oscillations. 
The oscillatory phenomena discussed here are para- 
metric excitation, synchronization, simultaneous oscil- 
lations and amplitude limitation by means of lamps. 
The reader should understand, however, that the meth- 
ods presented in this builetin can also be applied to a 
number of other oscillatory problems. 
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1. INTRODUCTION 


One of the most significant trends of the twentieth century is the 
departmentalization of science. The broad field of human experience 
has been divided into innumerable sections and subsections, each mas- 
tered only by a small group of specialists. One of the consequences of 
this division of science is that structural similarities between different 
fields are often overlooked. The recognition of these similarities can 
frequently be a great source of inspiration; a source closed to the 
specialist who fails to develop any interest outside his own field. 

A typical example of such a similarity is the occurence of oscillatory 
phenomena in such diversified fields as electronics, aerodynamics, eco- 
nomics, biology, etc. In recent years a uniform theory of oscillations 

applicable to all of these fields has been developed. However, the main 
applications of this theory at present are in electrical and aeronautical 
engineering. Mathematically, oscillatory phenomena lead to nonlinear 
differential equations. Because of this, the theory of oscillations is often 
called ‘“‘nonlinear mechanics.”’ 

Most important for the study of an oscillatory system are the steady- 
state oscillations. They may be stable or unstable, depending on whether 
the oscillator will or will not return to its original state after being 
subjected to a small disturbance (noise, ete.). An oscillator may have 

several possible steady-state oscillations. In this bulletin particular 
emphasis is placed on the nonsteady-state or transient oscillations. 
These terms, nonsteady-state and transient, are used synonomously. 
The study of these transient oscillations is often important for the 
complete understanding of the circuit behavior. For example, in a 
system having several possible steady-state oscillations, the transient 
study predicts which one of these will be reached from a given set of 
initial conditions. It also determines the manner in which oscillations 
will build up and describes the behavior of the system if a circuit para- 
meter is changed abruptly. Generally, a study of the transient oscilla- 
tions results in a more complete understanding of the oscillatory system. 

One important problem to be overcome by the investigator is that 
of representing the transient oscillations. Analytical representation is 
not practicable, since, even if an approximate formula could be found, 
it would be so complicated that it could not be interpreted easily. The 


7/ 
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most common graphical method of representation makes use of a plane 
with time as the abscissa and the dependent variable as the ordinate. 
This method is, however, limited to the cases in which the differential 
equation can be reduced to 


dx 


This reduction is possible for only a few oscillators. 
Autonomous systems with one degree of freedom can often be 
described by the differential equation: 


The corresponding oscillations can be represented in a plane with 2 as 
the abscissa and da/dt as the ordinate, the ‘‘phase-plane.’”’ Such a repre- 
sentation is shown in Chapter II. 

No general method for the graphical representation of oscillations 
in more complicated systems is available. If it can be assumed, however, 
that the oscillations are almost sinusoidal, then a large class of oscillators 
can be represented in some sort of a phase plane. In this new plane, an 
oscillation with a certain amplitude, phase angle, etc., will correspond 
to a single point. 


Il. AUTONOMOUS OSCILLATIONS WITH ONE 


DEGREE OF FREEDOM 


The study of the autonomous oscillators with one degree of freedom — 
takes a rather special place in nonlinear mechanics insofar as these 
oscillators can be treated conclusively and in all generality by topological 
methods. Equation 1.1, corresponding to autonomous oscillators with 
one degree of freedom, can be transformed into two simultaneous dif- 
ferential equations of the first order by the substitution dx/dt=y 


“Y= f(xy) y —9 (oy) 2 


dx g 


l| 


The time can be eliminated from these two equations by dividing 
one by the other: 


The solutions of this equation can be represented by curves in the a-y 
plane, which can be found by graphical construction (method of iso- 
clines). If the system is at rest, then « and y are constant in time; 
that is, 

dx _ dy 


Mie <hr 


and the differential, dy/dz, at such a point, is indeterminate (singular 


oint). 
4 The z-y plane may also contain a number of closed curves, or “limit 
cycles,” corresponding to periodic solutions of Eq. 1.1. The singular 
points and the limit cycles determine the structure of the x-y plane. 
The state of the oscillator is completely described at any one time by a 
point in the z-y plane. If the system is not disturbed, then this “‘repre- 
sentative” point will move along a curve defined by Eq. 2.2. 

Typical oscillatory systems that can be treated by this method are, 
for example, the ordinary triode oscillator, the Prony-brake and the 
multivibrator. It is possible but impractical to treat more complicated 
systems by this method because spaces with three or more dimensions 
are required to describe them. 


Ill. ALMOST SINUSOIDAL OSCILLATIONS 


A large group of oscillatory phenomena can be represented in a plane 
entirely different from that discussed in Chapter II provided that it 
can be assumed that the oscillations are almost sinusoidal. In electrical 
circuits, for example, oscillations are nearly sinusoidal if the quality 
factors (Q’s) of the resonant circuits are high. Then if it is assumed 
that the oscillations are purely sinusoidal, the circuits can be described 
completely by a few parameters; for example, by the amplitude of 
oscillation, the phase angle relative to an external sinusoidal voltage, 
etc. For steady-state oscillations, these parameters remain constant; 
for transient oscillations, they will change. A representation in a plane 
is possible if the system can be described by two parameters only, say 
X and Y. When one of these parameters is a phase angle, the plane 
lies on the surface of a cylinder so that the lines, ¢=0 and o¢=27, are 
identical. Many oscillatory phenomena can be described by two para- 
meters and can, hence, be represented in this plane. 

The rate of change of X and Y is determined by the state of the 
system and, hence, depends on X and Y only. + 


ax 
Ser ASN 
< = fo (X,Y) (3.1) 


Equations 3.1 can be obtained by various analytical and experimental 
methods. The method used in this bulletin is that of equivalent lineariza- 
tion. The mathematical formulation of this method has been described 
extensively in Buls. 395 and 400, and hence is not presented here. 

The time can be eliminated from Eqs. 3.1 by dividing one by the 
other. 


dX Aer) 
Was RG ee (3.2) 


The solutions of this differential equation can be represented as 
curves in the X-Y plane. If the system is not disturbed, then the — 
“representative point,” corresponding to an oscillation with parameters. 
X and Y, will move along one of these curves. For steady-state oscilla- 
tions, X and Y remain constant in time. 
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Dae av, 
if ee (33) 


_ This corresponds to a singular point in the X-Y plane. The symbols, 
X, and Yo, are used for the parameters of the steady-state oscillation. 
A singular point and the corresponding steady-state oscillation are 
stable if the oscillation and the corresponding point will return to their 
original position after a small disturbance. In order to discuss the 
behavior of the system in the neighborhood of the steady-state oscilla- 
tions, it is necessary to expand Eq. 3.1 around X_ and Yo. Letting 


6X = X — Xo 
Oe ue— 5 


then the first terms of the expansion are 


1AO25 beak eee 

7 Da Gn 
GG ee says fe 

Fe a ney atone (3.4) 


It is assumed that 5X and 6Y are small and therefore that the higher 
terms of this expansion can be neglected. Equation 3.4 is a set of two 
simultaneous linear differential equations which can be solved by the 
usual methods. The variables 6X and 6Y will approach zero from any 
initial conditions in the neighborhood of the steady-state oscillations 
provided that both the roots «: and x of Eq. 3.5 have negative real parts. 


Cygann oY kof 
a 
a oD) 


Necessary and sufficient conditions for this to exist are that: 


Of, Ofe 
VER py Oo 
Apri ies cate Of: 


LD Sea NCD 


20 (3.6) 


Like the limit cycles of the z-y plane described in Chapter H, the 
X-Y plane may contain closed curves which are called ‘‘drift curves.” 
Physically, these drift curves correspond to oscillations with periodically 
varying amplitude. 


IV. PARAMETRIC EXCITATION OF NONLINEAR SYSTEMS 


Oscillations in resonant circuits can be generated in several different 
ways; for example, by an external sinusoidal force or by a negative 
resistance. Another possibility is the excitation of the circuit by means 
_ of the periodic variation of a circuit parameter. This excitation is 

particularly strong if the capacitance or inductance of the circuit is 
varied with a frequency approximately twice that of the resonant 
frequency of the circuit. The amplitude of the oscillations so generated 
will build up to a value determined by the nonlinearities of the system. 

Assume for example that in the circuit of Fig. 4.1, the distance d 
between the plates of the capacitance C is varied according to 


d = d) +d; cos 2wt (4.1) 


where did) and w is close to the resonant frequency of the circuit. 
Then if the resistance R is sufficiently small, oscillations with a fre- 
quency w will build up. As is shown later, they may reach a stationary 
value if either R or L is a function of the current 7 (nonlinear resistance 
or inductance). 

In the first part of this chapter, it is assumed that R and L are 
constant; R= Ry, L =Lo. The capacitance C corresponding to Eq. 4.1 is: 


Co 


C= 
1+ y¥ cos 2wt 


where y =d,/d)<1. 

The circuit of Fig. 4.1 can be transformed to an “equivalent linear 
circuit” by the method of equivalent linearization described in Buls. 
395 and 400. The first step is to split the capacitance C into two parts, 
C, and Cy, as shown in Fig. 4.2. 


Cio = 


Co 
a 


wo? + ¥ cos 2uwt 


(4.2) 
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where w”=1/L Co is the resonant frequency of the system. It was 
assumed that ww ; therefore Ci is very large. 
For a large Cio and a small Ro, the current 7 is approximately 


t= TI cos (wt + ¢) 


The amplitude J and the phase ¢ will of course change, but over a few 
cycles they will remain practically constant. The voltage v across Ro 
and Cy) is (Fig. 4.2): 
p= : fiat + 1K 
cCn 
1 — w/w? — y/2 cos 2¢ 
wl 


fe 1( Re a 0,8 26 ) cos (wt + ¢) 


=I 


sin (wt + ¢) 


+1 0, sin (8wt + ¢) 
Due to the component of v at frequency 3 a small current at this 
frequency will flow through Ly) and C;. Since this component is very 
small, however, it can be neglected. As is shown, the components at 
frequency w will cause a change of J and ¢. 


Fig. 4.1. General Oscillatory Circuit Fig. 4.2. Oscillatory Circuit with Variable 
with Excitation by Variable Capacitance Capacitance Excitation. R and L Constant 


If all components of the current except those at frequency are 
neglected, then the circuit of Fig. 4.2 can be replaced by that of Fig. 


4.3 where: 


= x ; 
R.= Rot aC, sin 2¢ 


Co 
ots 1 — w?/a? — y/2 cos 2¢ oy) 


The circuit of Fig. 4.3 is linear and can be discussed by the method 
of linear analysis. 
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R Fig. 4.3. Equivalent Linearized 
Lo . Circuit of Fig. 4.2 


Te an 


The energy stored in the oscillator, LoJ?/2, will remain approxi- 
mately constant over one cycle. Over a longer period of time, it will 
change slowly since part of it is dissipated in the resistance R,. 

d (IyI?/2) [? 


dt 7g. te 


or 


CLF ay Ben cA ) 
dt = 2L, (fa + 52 sin 2 (4.4) 


Similarly, the phase angle ¢ will change slowly. 


dp — 1 
@ as J 7 CC. 
. Cy > oF 

or 

do 1 G 

di e520 CA ph 

= Aw — f © cos 26 (4.5) 

where 


Aw = Dy <=> oF 
The time can be eliminated from these two equations by dividing one 


by the other: 
dp _ Aw — (v/4)w cos 2¢ 


~ (E/2La)[ Ro +32 sin 26 | 
0 


2 — cos 2p 
—I(p + sin 2¢) 


(4.6) 


where 


p= 


2 is proportional to the “detuning” of the resonant circuit. 


PR ah wn 
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The trajectories in the J-¢ plane corresponding to Eq. 4.6 permit 
graphical representation of the transient solutions. They can be ob- 
tained either by the method of isoclines or by direct integration. Typical 
groups of trajectories are shown in Fig. 4.4. 

Steady-state oscillations exist if dl/dt=dd/dt=0. This corresponds 
to singular points in the [-¢ plane. The stability of these singular points 
and of the corresponding steady-state oscillations can be determined 
from the behavior of the trajectories in their neighborhood. It can be 
seen that they are unstable, because after a small disturbance of ¢ the 
amplitude J will increase or decrease indefinitely. 

In order to obtain stable steady-state oscillations in the circuit of 
Fig. 4.1, it is necessary that the circuit contain a nonlinear inductance 
or a nonlinear resistance. 


g g 
7 277 
37 3m 
a 2 
Zz 7 
qn q 
4 2 
Z / 
(a) z=0, p=O 0 (b) -2=2,.p =0 
g g 
| ae 
51 
7H clone 
uh 
Z Re 
1T 
A é 


aria 


O(c) z-£, p=F# (Small Resistance) O (@) z=, p=2 (Large Resistance) 


Ds 


Fig. 4.4. Trajectories 
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Consider first the case where the inductance L of Fig. 4.1 is a fune- 
tion of the current 7 passing through it: 


L=btl te) “<1 


The oscillatory curcuit corresponding to Fig. 4.1 with a nonlinear in- 
ductance is shown in Fig. 4.5. It is again assumed that 


zt =I cos (wt + ¢) 


The voltage v is then: 


d : : 1 : 
v= a (Loe*) + Rot seas tdt 
3 1 — w/w? — y/2 cos 2 : 
= 1 |—elol*a + let AL? 008 PFI in (wt + 6) 


rad (Re t SHEP cos 26 ) cos (wt + ¢) 


+ terms at frequency 3w. 


The terms at the frequency 3w are again neglected. The equivalent 
conductance and capacitance are: 
= seks ae 
Re. = Ro + aC, sin 2¢ 


Cr Gs (4.7) 


1 — w/w? cos 26 — Sop 


The differential equation corresponding to Eq. 4.6 is: 
do Aw —(y/4)w cos 26 — 3% eP2w 
dl — (1/2Lp)[Ro + (y/2wC>)sin 2¢] 


_ 1% — cos 2¢ — al? 
I (p + sin 2¢) 


where 
ep aS) Aw 
woy 
2RowC 
Vf 
eonrteats EPES: 
24 
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Fig. 4.5. Oscillatory Circuit with 
Variable Capacitance Excitation. 
Inductance Nonlinear 


For a steady-state oscillation it is necessary that 


GE AG 

GE Reread Pe 
or or! (0), 

z — al? — cos 2¢ = 0 


p+sin 2¢ = 0 (4.8) 
The symbols J) and ¢o are used for the values of J and ¢ that satisfy 


Eqs. 4.8. For each set of z, a, p there are at most two independent 
solutions: 


ae: ae 2 
I? = = 2 (cos 2¢ > 0) 
and 
an ee (cos 2¢ < 0) 


Only one of these solutions is stable, however. 
The characteristic equation corresponding to Eqs. 4.4 and 4.5 is 


x To OR. cua Io Ok, 
wCo = 1 ) ae (4 -) es 
Foe ol, ‘ 


Necessary and sufficient conditions for stability are that 


I, OR. Go (4 -) 
we <0 
= al 2 


ar jee fot ina -)>0 
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UMhichhdl Sato, 
Ege or Fig. 4.6. Region of Stable Oscillations 
table 
f 0) 
Oscillations oe 
+/ 
SS 
P 


RS) 


Ny 


Fig. 4.7. Trajectories for p =0, z= 2 (Large Detuning) 


3 To i eae ee 


NIy 


Fig. 4.8. Trajectories for p = "%), Z= 2 (Large Detuning) 
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or, in terms of ¢ that 
sin 2¢ < 0 
a cos 2a < 0 


The first of these inequalities is satisfied if p>0O (Eq. 4.8). The 
second inequality shows that, of the two solutions for Jo, only one is 
stable. For a>0 the solution with cos 2¢<0 is stable; for a<0 the 
solution with cos 2¢>0 is stable. The regions of stable oscillations can 
be represented in the p-z plane (Fig. 4.6). Additional boundaries of the 
region correspond to: J? > 0 and p = —sin 2¢ < 1. The region does 
not depend on the magnitude of a as long as a remains finite. If p is 
small, then the inductance may arbitrarily be increased or decreased 
depending on the sign of a. This result, however, is due to the many 
simplifications made for this calculation; for example, with regard to 
the dependence of the inductance, L, on the current, 7. In any physical 
system, stable oscillations will occur only within a finite range of varia- 
tion of the inductance. Figures 4.7 and 4.8 show two sets of trajectories 
corresponding to systems with nonlinear inductance. For Fig. 4.7 the 
system is assumed to contain no damping (p=0) and large detuning 
(¢=2). The singularity is of the center-type and all trajectories are 
also drift curves. The case p=0 cannot occur in a physical system, 
however, since a small amount of damping is necessarily contained in 
any system. Figure 4.8 corresponds to small damping (p = 4) and large 
detuning (z=2). The stable singularity is a focal point which, however, 
cannot be reached by the system from rest (J =0) unless the detuning 
is considerably smaller. This corresponds to hard self-excitation of the 
ordinary oscillator. 

Stable steady-state oscillations may also exist if the circuit contains 
a resistance that is a function of the current 7 passing through it. 
(See Fig. 4.9.) 

R= Rj (1. + uv’) 


Ry and » may be positive or negative. If Ry <0, u>0, and y=0 (C is 
constant in time), then the circuit of Fig. 4.1 will oscillate with an 
amplitude 
| [v= ate (4.10) 
which result is used later. 

It can be shown (Fig. 4.10) that the equivalent capacitance and 
resistance are 

Co 

1 — w?/wo? — (y/2 cos)2¢ 


C. = 


R. 


Ro(1 cia oe sin 26 
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R is Re 
C Cio Gq Ce 
[Eales Seat g Laat 
Fig. 4.9. Oscillatory Circuit with Variable Fig. 4.10. Equivalent Linearized 
Capacitance Excitation. Resistance Nonlinear Circuit of Fig. 4.9 
and that the differential equation corresponding to Eq. 4.6 is 
do sis a Ore (4.11) 
al I(o(1+ 34 uF) + sin 24) 
where 
A 
Hes 
wy 
2RowC 
p = 
ti 
The conditions for steady state oscillations are 
z— cos 2¢ = 0 
p(1 + 34ul?) + sin 26 = 0 (4.12) 


and the conditions for stability 


sin 26 < 0 
pu > 0 : 
Figure 4.11 shows the region in the p-z plane where stable oscilla- 
tions are possible. Of these stable oscillations two different cases must 
be distinguished — p <0 and p>0. 


Fig. 4.11. Region of Stable Oscillations 


Bul. 421. OSCILLATIONS IN NONLINEAR SYSTEMS 21 
(a) Oe pel u>O0 


The amplitude of oscillation J can be found from Eqs. 4.12. The 
_ dependence of the amplitude on detuning is shown in Fig. 4.12. Outside 
the “region of synchronization” the amplitude of oscillation is zero. A 


group of trajectories for z=14, p=, w=1 is shown in Fig. 4.13. 
(b) p <0 <0 


The dependence of the amplitude on the detuning for this case is 
shown in Fig. 4.126. Inside the region of synchronization there is little 
difference between this and case a, but outside the region the amplitude 
does not drop to zero. This is due to the fact that even without para- 
metric excitation the circuit will oscillate. 


AL 


(a) p>O (6) pO 


Fig. 4.13. Trajectories for Circuit Nonlinear Resistance, z= 1/2, P = 13, R= 1 
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Fig. 4.15. Oscillations Corresponding to Drift Curve Shown in Fig. 4.14 


Figure 4.14 shows the trajectories for z=1.05. Because this is out- 
side the “region of synchronization,” no singular points exist. There 
exists, however, one (and only one) drift curve shown dotted in Fig. 
4.13. Figure 4.15 shows the oscillation corresponding to that curve. 


This figure shows a good agreement with some experimental results 
published by W. L. Barrow.* 


*Parenthesized superscripts refer to correspondingly numbered entries in the Bibliography. 


_ V. SYNCHRONIZATION 


In this chapter the oscillatory circuit of Fig. 5.1 is discussed. This 
circuit is the equivalent of an ordinary plate-tuned oscillator with an 
external sinusoidal synchronizing voltage v,. When »v, is zero, the circuit 
will oscillate with a frequency, 
re eee Da Ie 
a eLGe 
The frequency of oscillation w for v.40 is not necessarily w. As the 
frequency of the external voltage w, is varied and becomes approxi- 
mately equal to wo, w will suddenly change from w to w,. This phenomena 
is called locking, synchronization, or entrainment of frequency. It occurs 
not only if w,2w) but also if w,/w.=p/q, where p and q are “small 
integers.” This phenomenon is called subharmonic synchronization and 
has been discussed in detail in Bul. 400. 

If it is assumed that v, = V; cos wst, then the voltage v; between points 
A and B is approximately sinusoidal, particularly if the losses in the 
resonant circuits are small 


v1 & V cos (wi + 6) 


®o 


It was shown in Bul. 400 that the circuit of Fig. 5.1 can then be 
reduced to that of Fig. 5.2 in which the non-linear element is replaced 
by an “equivalent network.” The differential equations for V and ¢ 
for this circuit are 


nV 


nrging —9¢@ |G + Gel (5.1) 
do ak Cc 
SS — = w)— 5.2 
y= Ae Ty 3 (5.2) 
where Aw-= wo — LESS 
Pp 


These equations determine the transient behavior of the oscillator. 
Dividing Eq. 5.2 by Eq. 5.1 eliminates the time 


Vee OG AAG.) 
This. equation permits representation in a plane with V and ¢ as 
coordinates, the V-¢ plane. Steady-state oscillations (dV /dt =do/dt =0) 
correspond to singular points. Besides singular points, the V-¢ plane 


dod 2AwC a aCe (5.3) 


23 
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4=Veos(wt to) 
=Kcosat 
peter ct, 


Fig. 5.2. Equivalent Linearized Circuit of Fig. 5.1 


may also contain drift curves. They occur particularly if an oscillator 
is close to synchronization.” 

As an example, the special case w =3w, is presented here. The fre- 
quency of the oscillator is three times that of the external voltage 
(frequency-multiplication). It is assumed that the current 7 of Figh}a 
can be expressed by a third degree polynomial of the voltage v 


t= av + Bo? + yo? 


The equivalent conductance and capacitance for w=3w, and a third 
degree polynomial were calculated in Bul. 400. They are: 


ESP a SNe rR ee Rela a 
ex tesanisiate tivaee | Lis) ster py cos ¢ 
Re aoe ae 
Cis Ieoy VSm o (5.4) 


Using these two expressions, Eq. 5.3 may be written as 


dé __—-2Vi2 +(V8/V) sin o | 
RES AG ES (5.5) 


where 


Introducing the dimensionless parameters X=V/V, and X-/Vo, Eq. 
5.5 can be simplified to 


io 22 Xz+ X3 sin ¢ 
dX NBR (X= 1) eh eps) (5.6) 


The solutions of this differential equation can be represented as 
trajectories in the X-¢ plane. The singular points in this plane are 
stable if they fall inside the region of stability in Fig. 5.3. Figure 5.4 _ 
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shows the family of trajectories in the X-¢ plane corresponding to 
a stable singular point (point A). This point is a focus. Figure 5.5 
shows the trajectories in the X-¢ plane for an unstable singular point 
(point B). This figure also shows a drift curve corresponding to “near 


synchronization.” 
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Fig. 5.5. Trajectories with Unstable Singular Point 


VI. SIMULTANEOUS OSCILLATIONS 


One of the most controversial and also one of the more interesting 
phases of nonlinear mechanics is that of simultaneous oscillations in 
circuits with two degrees of freedom. In recent years a number of 
papers dealing with this problem have been written and it has been 
proven both theoretically and experimentally that simultaneous oscilla- 
tions are possible.“”"*” 

For simultaneous oscillations, the voltage between points A and B 
of Fig. 6.1 is, approximately, 


V1 = Vi cos (wit + 1) 
and the voltage between B and C is 
vg = V2 cos (wot + 2) (6.1) 


where w,4o . A further distinction between w; and a gives rise to the 
necessity for distinguishing between synchronous and asynchronous 
simultaneous oscillations. For synchronous oscillations the frequencies 
@, and w: are related by po1=qwe, where p and gq are small integers, 
and comes about as a result of an effect of mutual synchronization. 
When the sum (p+q) is large, the synchronization becomes so weak 
that it can no longer counteract the continuous disturbances caused by 
line fluctuations, noise, etc. For asynchronous oscillations no relation 
pw =qw. is preserved over an appreciable length of time. Mathemati- 
cally it is assumed in this case that the ratio w/w. is constant but irra- 
tional. Asynchronous simultaneous oscillations were discussed _ briefly 
in Bul. 395. 

The oscillations in a circuit with two degrees of freedom can be 
represented in a four-dimensional phase space with (v1, dv: /dt, v2, dv, /dt) 
as coordinates. For a successful treatment, however, it is imperative 
that the oscillations of any system be represented in a plane or at most 
in a three-dimensional space. The four dimensions of the present case 
ean be reduced to three if it can be assumed that »: and t are sinusoidal. 
The oscillations can then be represented in a space having as coordinates 
Vi, V2 and, for synchronous oscillations, the angle ¢=p¢i1—q¢2 which 
is defined as the phase angle between v and v. An oscillation which 
has amplitudes V1, V2 and phase angle ¢ corresponds then to a point 


in this space. 
27 
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oC 
Fig. 6.1. Oscillator with Two Degrees of Freedom 


The voltage v of Fig. 6.1 is 


v =v, +» = Vicos (wt Se $1) a V2 cos (wot —- de) 
This voltage applied to the nonlinear element will produce a current 
0 = In cos (wnt + $1) ot Tie sin (wt + $1) 


+ I21 COS (wet + 2) + Io2 sin (wet + dy») 
+ terms at frequencies other than «, or w». (6.2) 


As shown in Bul. 395, Chapter I, the circuit of Fig. 6.1 can then be 
replaced by the two circuits of Fig. 6.2 where 


Goats Gem 8 
: thee (6.3) 

mee Ce AE eco! 

Ge si V2 ) Core a eis 


where the equivalent conductances and capacitances are functions of 
Vi, Vo, and ¢. 

If p and q are small integers, then the expression o=p¢d1— 2 is 
determined independent of a shift of the time axis. However, if either 
p or q is irrational, ¢ can still be adjusted arbitrarily close to any value 


Ge trequency uy 


Ge Trequency w, 


> 


Fig. 6.2. Equivalent Circuit of Fig. 6.1 
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by a suitable shift of the time axis. Since this shift is a purely mathe- 
matical operation and cannot have any influence on the oscillator, it 
is clear that the equivalent impedances for asynchronous oscillations 
- cannot be functions of ¢. 

The transient and steady-state behavior of the system can be 
calculated from the circuits of Fig. 6.2. The conductances of Fig. 6.2 
will damp the oscillations and change the amplitudes. The capacitances 
C; and Cy, will change the frequencies, w; and we, and cause a shift of 
the phase angle, ¢. 

The corresponding differential equations are for asynchronous oscil- 
lations 


Over Vi 
ia IG 
CN Se eM 
Se ee oC (Go + Gee) (6.5) 
and for synchronous oscillations 
OV peer 
OVet 2 Ve 
dt 7 IC; (Ge Sia Ge) (6.7) 
ag Pon Ce 
a Aw ante (6.8) 
where 
Aw = P d 
“/ IhC, / L2Cs 
C 
Ge — Cu Fa G, Ce 


By dividing Eqs. 6.5 it is possible to eliminate the time between them, 
giving 
dV, — ViC2 (Git Gre) 


= 6.9 
dV» VoCi (Go + Gre) 66.2) 


This differential equation can be solved by the methods of isoclines 
in a plane with Vi and V2 as coordinates. A point in this plane corre- 
sponds to an oscillation with the amplitudes, Vi and V2. If the system 
is not disturbed, then the representative point will move along a tra- 
jectory. Figures 6.3 and 6.4 show typical groups of trajectories. As 
shown in Bul. 395, stable simultaneous steady-state oscillations are not 
possible if the current and voltage across the nonlinear element are 


related by 
i = av + Bo? + yw (6.10) 
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Fig. 6.4. Trajectories for Asynchronous Oscillator and i= av + Bv? + yv? + bv! + ev? 
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The trajectories shown in Fig. 6.3 correspond to those obtained from 
such a polynomial. The point A shown in this figure is one that would 
ordinarily correspond to simultaneous oscillations, but it is a saddle 
point and is therefore unstable. On the other hand, if the current- 
voltage relationship is a fifth-degree polynomial and if certain other 
conditions are satisfied, then stable simultaneous oscillations are possi- 
ble. Such a condition is represented by point B in Fig. 6.4. This figure 
also shows that as a rule asychronous simultaneous oscillations are not 
self-starting. 

In general the representation of the synchronous oscillations in the 
V,-Vo-@ space is not practicable. For certain conditions, however, 
Eq. 6.8 can be reduced to 


£6 =f (Vi,Vs) sin @ (6.11) 


where f(Vi, V2) >0 for Vi>0, V2>0. All trajectories will then approach 
the ¢=7 plane and a trajectory starting in this plane will remain in it. 
Figure 6.5 shows a typical group of trajectories in the ¢=7 plane. 
The simultaneous oscillations corresponding to point A in this figure 
are stable and self-starting. 


Fig. 6.5. Trajectories for Synchronous Oscillations (#1 = 32) 


Vil. THE LIMITATION OF THE AMPLITUDE OF 
OSCILLATIONS BY LAMPS 


The amplitude of oscillations in an oscillator containing only linear 
elements will either decrease or increase indefinitely. It can reach a 
stationary value other than zero only if the oscillator contains a non- 
linearity of some sort, in most cases an electron tube. The electron tube 
has a double function: it delivers a-e power to the passive elements 
and limits the amplitude of oscillation. 

In the oscillator discussed here, these functions have been separated. 
The electron tube delivers the a-c power, but the amplitude is limited 
by a lamp inserted in the resonant circuit. The temperature of this 
lamp depends on the power dissipated in it and, hence, on the amplitude 
of oscillation. The resistance of the lamp varies as a function of fila- 
ment temperature. 

For small amplitudes of oscillation more power will be generated in 
the electron tube than is dissipated in the lamp, and as a consequence 
the amplitude will increase. Due to this increase, the temperature 
of the lamp, and also its resistance, will become greater. At a certain am- 
plitude the resistance of the lamp will be so great that the power gener- 
ated in the electron tube will be equal to that dissipated in the lamp, 
and steady-state conditions for amplitude are achieved.” 


Fig. 7.1. Oscillator 


An idealized oscillator using a lamp to limit the amplitude of oscilla- 
tion is shown in Fig. 7.1. Actual oscillation using lamps to limit their 
amplitudes are usually more complex.” 


The variation in resistance of a conductor with the temperature is 


Rr, = Rr, (1 + aT) (7.1) 
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where @ is a constant depending on the material; Rr, the resistance at 
the ambient temperature 7); and T;=T—T) the difference between 
the temperature of the conductor and the ambient temperature. 

It would be difficult to treat the circuit of Fig. 7.1 in full generality; 
therefore a number of simplifying assumptions must be made. It is 
assumed that no grid current flows and that only the linear portion of 
the triode characteristic is used. The plate-current is then a function 
of the plate voltage only: 


: il 
tp = Jm\ Vg + — v») 
be 


1 
a om (4 =) Vp = —GwWp Ceo) 


. where 


iota: 
is a constant with the dimensions of a conductance. Equation 7.2 shows 
that the triode of Fig. 7.1 can be replaced by the equivalent negative 
conductance —G». (Fig. 7.2) The lamp has been replaced in this figure 
by a resistance fr,. 

The circuit of Fig. 7.2 can be further simplified. Since the frequency 
of oscillation w differs by only a very small amount from 1/,/ LC, 
the conductance —G) can be replaced by a resistance — fp in series 
with the inductance L (Fig. 7.3). 


Ro = Go 


The total resistance of this resonant circuit is 


R= —Rk) + Rr, 
= —It + Rr, (1 + a1) 
= —R,(1 — BT:) (7.3) 
where 
Ri =k) — Rr 
Rr, 
Caan 


It is assumed that the quality factor Q of the circuit is high. The 
current 7 can then be approximated over a few cycles by 


t =I cos at (7.4) 
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Fig. 7.2. Oscillator with Triode Replaced by Negative Conductance 


The amplitude of current J will vary slowly and, for a short period of 
time, it can be assumed constant. 

The total energy stored in the circuit of Fig. 7.3 is LJ*/2. During 
one cycle, it changes twice from electromagnetic to electrostatic energy 
and back again, but the total amount remains approximately constant. 
A slow change is caused, however, since some of the energy is dissipated 
(or generated) in the resistance R. 


Le di? Maer, 
Os df eae 


I2.(1 — pT;) (7.5) 


A similar differential equation can be found for the thermal energy 
mcT, of the lamp 

dT, 2 i? 

eee —-KT™%+ 5 

where m is the mass of the heat element, c is the specific heat of the 

element, and K is the constant of thermal conductivity. For steady- 

state oscillations it is necessary that both J and 7, remain constant; i.e., 


Rr, (7.6) 


Pie nae Tae (7.7) 
These oscillations may be stable or unstable depending on whether they 
will or will not resume their original amplitude after being subjected to 
a small disturbance. This can be determined by calculating and examin- 
ing the conditions for stable steady-state oscillations from Eqs. 7.5 and 
7.6. Such an examination shows that these oscillations are stable. 


Fig. 7.3. Equivalent Circuit of Fig. 7.2 
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The behavior of the nonsteady-state oscillations can be determined 
also from Eqs. 7.5 and 7.6. Dividing one of these equations by the 
other eliminates time as a variable and makes possible the representation 
of the transient oscillations in a plane with J? and 7; as the coordinates 


aly EL “OK [Ro sae hy (1 a B ss T,)| 


dP ~ Rit Pd — 6st, Wess 


where 7 = K/mc is the thermal time constant of the lamp. 

This differential equation cannot be solved in closed form. An ap- 
proximate solution, however, can be found by the method of isoclines. 
In Fig. 7.4, a typical solution is represented by a group of trajectories 
in a plane with the amplitude of oscillation as the abscissae and the 
temperature as the ordinates, the J°-T; plane. 

At a given time the oscillator will have an amplitude, J, and the 
lamp a temperature, 7). The point in the J?-T, plane defined by these 
two values will be called the “representative point” since it represents 
the state of the oscillator. As temperature and amplitude vary, the 
representative point will move along a line or “trajectory” in the 
[?-T; plane. For example, if the oscillator starts from rest (point 0 
in Fig. 7.4), it will move toward the point B along the dotted line. 
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Temperature 


0 
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Fig. 7.4. Trajectories Describing the Transient Behavior of the Oscillator 
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The point B corresponds to steady-state oscillations with amplitude, 
I,, and temperature, T,. 

The steady-state oscillations corresponding to point B will be dis- 
turbed constantly because of noise, variations in circuit parameters, 
and change of the ambient temperature; but it can be seen that the 
representative point, after having been displaced from B, will always 
return to its original position. Consequently the steady-state oscillations 
corresponding to point B are stable. 

By a similar consideration it can be seen that point A is unstable 
since the representative point will never return to this point after it 
has been given a small displacement. 

By assuming the current 7 to be 


2 =I cos at 


it was possible to describe the oscillator by the equations 7.5 and 7.6. 
If no assumption regarding the current 7 is made, however, it is neces- 
sary to use three differential equations of the first order. For example, 
the circuit of Fig. 7.3 leads to 


di 
wae 
d 1 1. 
dT 7 1 
~ = ar + —— (Ry — Ri (1 — BT) 2 (7.9) 


dt fh mc 


where 7 is the current through the resonant circuit. 

Similar to the [?-T; plane it is possible to set up a 7-y-T; space. 
The solutions of Eqs. 7.9 could then be represented by a group of 
trajectories in this space. 


Vill. CONCLUSIONS 


Several oscillatory systems have been discussed in this bulletin with 
special emphasis on nonsteady-state or “transient” oscillations. For 
autonomous systems with one degree of freedom and any degree of 
nonlinearity, a complete representation of the transient oscillations is 
possible in the phase plane. For more complicated systems, this method 
of representation is either not possible or impracticable since it leads 
to spaces with four or more dimensions. 

A representation in a plane (which is no longer the phase plane) is, 
however, still possible in some important cases if it can be assumed 
that the oscillations are almost sinusoidal. As has been shown, the 
transient oscillations can then often be described by two first order 
differential equations: 


dx 
ergs Jt aah, 
ay: 
Serre (X5¥) (8.1) 


where X and Y may be an amplitude, phase angle, temperature, etc. 
These equations do not give a complete description of the oscillations, 
however, since they neglect the higher harmonies. This does not lead 
to appreciable errors if the losses in the resonant circuits are small 
(high Q) since then the amplitudes of the higher harmonics are also 
small. Even for low Q circuits the results calculated from Eq. 8.1 are 
confirmed surprisingly well by experiments. 
The general procedure is to eliminate from Eqs. 8.1 by division: 


dX _ filX,Y) 
dY — fo (X,Y) 


The solutions of this differential equation can then be represented as a 
group of trajectories in the X-Y plane, in which an oscillation with 
parameters X and Y corresponds to a point lying on one of the tra- 
jectories in this plane. If the system is not disturbed, then this point 
will move along the trajectories defined by Eq. 8.2. Steady-state oscil- 
lations correspond to singular points of this equation. These oscillations 


(8.2) 
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are stable if after being subjected to a small disturbance the parameters 
X and Y will return to this point. Besides the singular points, the X-Y 
plane may also contain closed trajectories or drift curves. 

Four different oscillatory phenomena were treated by this method: 
parametric excitation, synchronization, simultaneous oscillations and 
amplitude limitation by means of lamps. The method can be applied 
to a number of other problems, such as oscillations in conservative 
systems with several degrees of freedom, self-modulating oscillators, and 
the motor-boating effect occurring in audio oscillators. 

The results obtained in Chapter V and VI were confirmed experi- 
mentally. The degree of correspondence between theory and experiment 
was good. 
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